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Abstract 



CN \ We study the asymptotic stability of de Sitter spacetime with respect to 

— ^ | non-linear perturbations, by considering second order perturbations of a flat 

£> ■ Robertson- Walker universe with dust and a positive cosmological constant. 

~D ■ Using the synchronous comoving gauge we find that, as in the case of linear 

perturbations, the non-linear perturbations also tend to constants, asymp- 
totically in time. Analysing curvature and other spacetime invariants we 
show, however, that these quantities asymptotically tend to their de Sitter 
values, thus demonstrating that the geometry is indeed locally asymptoti- 
cally de Sitter, despite the fact that matter inhomogeneities tend to constants 
in time. Our results support the inflationary picture of frozen amplitude mat- 
(2JT), ter perturbations that are stretched outside the horizon, and demonstrate the 

validity of the cosmic no-hair conjecture in the nonlinear inhomogeneous set- 
tings considered here. 
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The cosmic no-hair conjecture |3J states - roughly speaking - that 'all expanding uni- 
verse models with a positive cosmological constant asymptotically approach the de Sitter 
solution' 0|. Previous studies of this conjecture fall into two main groups: those involving 
special models, e.g. containing Killing vectors or special Petrov types, and those using a 
linear perturbative analysis. Among the first group is the work of Wald 0, who proved 
the conjecture for ever expanding homogeneous models, as well as studies of inhomogeneous 
models J3]. Among the second group are studies of the stability of de Sitter with respect to 
linear tensor (gravitational wave) perturbations [|J, as well as investigations involving both 
linear scalar and tensor modes 0,[J, where the scalar modes arise in presence of matter. 
The authors of [|5],|6| studied the dynamics of linear perturbations in Friedmann-Lemaitre- 
Robert son- Walker (FLRW) models with dust and a positive cosmological constant A, by 
employing a synchronous comoving gauge, and showed that both scalar and tensor modes 
tend to constants in time, while the vector modes were shown to decay ||. This establishes 
the stability of de Sitter in these linear perturbative settings, but not its asymptotic stability 
||. It was also shown in ||, however, that there exists a coordinate transformation (valid 
within the horizon, see below) in terms of which the asymptotic metric becomes locally 
de Sitter in static coordinates ]7). This at least partly demonstrates that one in fact has a 
local asymptotic stability, in the sense that for a freely falling observer in such a frame the 
spacetime approaches de Sitter in time, exponentially fast. 

A number of studies have addressed this conjecture from a broader perspective (see e.g. || 
and Refs. therein), and in particular linear perturbations of de Sitter have been studied 
extensively in the context of the inflationary scenario (see e.g. ||[10| and Refs. therein). 



The latter has been one of the main motivations for the interest in the cosmic no-hair 
conjecture. Indeed inflation was devised as a mechanism to account for the observed large- 
scale homogeneity and isotropy of the universe, and while it is natural to assume a FLRW 
model to work out observational predictions of inflationary models such as perturbation 
spectra, it is nevertheless important to address the fundamental question of whether locally 
a FLRW geometry can arise from an open set of generic initial data, whose memory is lost 
through an inflationary phase. 

The aims of this Letter are twofold: to extend the perturbative analysis of the stability 
of de Sitter to second order and to provide an invariant characterization of the asymptotic 
results of both first and second order perturbations. 

Firstly, in order to make contact with earlier results, we study the asymptotic evolution 
of the first and second order perturbations of a flat FLRW model with dust plus A. Asymp- 
totically, this can be viewed as a de Sitter universe with dust perturbations. Secondly, to 
throw more light on the asymptotic state of the perturbed universe model, we employ cur- 
vature and other spacetime invariants. This Letter contains a short account of our results; 
a more detailed analysis will be presented elsewhere flll |. 



We proceed by using the formalism of [|T2],|T3|, extending it to include A. In particular we 
work in the comoving synchronous gauge and assume vanishing vorticity. There is no loss 
of generality in the latter choice, as vorticity decays during expansion. The line element of 
the perturbed spacetime can be written as 

ds 2 = a 2 (r)(-dT 2 + - 1a ^dx a dx (i ) , (1) 

where a is the scale factor, r is the conformal time given by dr = dt/a(t) and, in usual 



notation 
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7^ = ^ + 7^ + ^- (2) 



The full equations to be considered are: those for the evolution of the expansion tensor 



a (3 



O'a/3 + Q@al3 + R *a/3 = ^P^afi + A<5 Q/3 , (3) 



its trace, the Raychaudhuri equation 



& + o a(3 o a p + l - P = a , (4) 

the continuity equation 

p' + 0p = 0, (5) 

and the energy constraint equation 

e 2 -e^e Ql3 + R* = 2(p + A), (6) 

where p is the matter-density, Op = u"p = ^■y aSr y'sf3> R* * s the 3-Ricci scalar and the prime 
denotes the derivative with respect to r. Since we are interested in the asymptotic behavior 
of the perturbations, we shall assume that the background evolution has already reached the 
asymptotic regime, with the A term dominating in the Friedmann equation. The Friedmann 
equation can in this case be solved to give 

<r) = ~, H=^, (7) 

where the conformal time r is negative and tends to zero as a goes to infinity. 

We start by considering the first order perturbations since they act as source terms in 
the second order perturbation equations. In the linear theory the scalar and tensor modes 
are independent and for an irrotational spacetime one can set the vector modes equal to 
zero. The evolution equation for the scalar perturbations is obtained by combining (H), © 
and @ and can be written in terms of the density contrast 5 = (p — pb)/p in the form 

S" + £$' _ la?p§ = o , (8) 

where p = po/a 3 is the background density and po is a constant. Note that the last term 
of this equation is asymptotically negligible and therefore the asymptotic solution of (||) 
has two modes, one constant and the other decaying. Expanding the full solution of this 
equation around r = one obtains up to second order in r 

*(r,x) - CMx) tej£ + l^r 2 ] + C 2 (x)r 2 , (9) 



where A 2 = hpoH; C\ and Ci are arbitrary spatial functions related to the initial data. The 



scalar metric perturbations are related to S |12] and also tend to constants in time |l I 



The tensor perturbation K a p is the transverse traceless part of ^A and is gauge- invariant. 
The asymptotic form of the equation for 7r a/ g is similarly obtained from the linearization of 
the evolution equations, and for each Fourier mode k is given by 



7T. 



2 , 



a ,3 ~ ~<0 + frXtf = . (10) 

r 

This equation can be solved to give a solution for each mode k 

^a/3 = a a/3(sin kr — kr cos kr) + 6 a/3 (cos kr + kr sin kr) , (11) 

where a a p = a a/ 3(k) and b a/3 = bapik). A formal solution ^^(x, r) can then be obtained in 
terms of a Fourier integral, which can asymptotically be approximated by 

(12) 

Thus both the first order scalar and tensor perturbations @ and ( |i~l~D asymptotically ap- 
proach constants in time. 

The analysis so far confirms known results ||-f| and sets the stage for proceeding to the 
second order. Before doing so, it is however worthwhile taking a brief look at some first 




order gauge-invariant quantities. The Electric Weyl tensor satisfies pi 



a®V b E h a =\pV a & 2^<5> H = pe m , (13) 

where the equivalence with the equation on the right hand side (Eq. (4.3) in Bardeen [[0J) is 
in Fourier space; T> a is the comoving fractional density gradient of density, a covariant gauge- 
invariant measure of inhomogeneity |16|| , and e m is the equivalent Bardeen variable [|i4" 



Now both these quantities evolve exactly as the density perturbation in the comoving gauge 



above, and thus tend to constant values in time in the present context. The simple Eq. (13 



is the relativistic perturbation equivalent of the Poisson equation, and Bardeen's potential 



$# is the gauge-invariant analogue of the Newton's potential |14| . This equation carries a 
great deal of useful information. Firstly it shows that in the absence of matter, the Electric 
Weyl tensor is transverse, so that there are no scalar modes in this case, $# = 0, but only 
gravitational waves (cf. [|17]]. Secondly it shows that $# decays as ~ or 1 when the density 
perturbation is frozen, with a corresponding decay in (the scalar part of) E h a . Finally $# 
directly corresponds to the metric perturbation variable in the longitudinal gauge || , which 
therefore also decays (cf. fl2|1 ). 

Concerning the second order perturbations, the spatial metric in (|TJ) can be written as 

7 g = -20^ a/3 + X g , (14) 

where <fi^ and x a p denote respectively the trace and trace-free parts. The usual procedure 
is then to substitute ( |T4D and (0) in the field equations, using a similar procedure as in the 



linear case. Keeping the lowest order terms of @ and (|P2|) on the right hand side of the 
Raychaudhuri equation we obtain 



(2) " - V' + \^P^ {2) = F (*)t + 0(r 2 ) , (15) 

where the source term F is a spatial function (given in [|ll|]) that depends quadratically on 
first order perturbations, and in particular from the four initial data quantities Ci,C2,a a p 
and b a p. The solutions to fll5|) have the form 

(2) = ^ + rC 3 (x) J f (3^) + rC 4 (x)F § (^) , (16) 

with C3 and C4 are arbitrary spatial functions. The asymptotic behavior of this solution 
can be obtained by expanding around r = to give, up to second order in r 

(2) -^ + ^3(x)^j+C 4 (x)r 2 . (17) 

Now, assuming standard second order initial conditions [|12] 0^(ro,x) = and (f)^ 2 ' (to,x) = 
we find 



~ F I A2 n b£) ~ F / A2 

;f - B 2 )tZ - Bl ' UW 1 + (f - l)r 2 



C 3 (x) = -= '— , C 4 (x) = -^ — , (18) 



where i?i = r(|)/(y437r) and B 2 = 3eA3/(6r(|)). Therefore, 0^ 2 - ) asymptotically approaches 
a constant value in time which depends on r , C±, C 2 , a a p and 60/3. Furthermore, from the 
momentum and energy constraints, we obtain up to second order in r 

x£V, x) « g,(x) , X <$*( T , x) « G(x) , (19) 

where Qp and G are spatial functions given in [D]] . Finally, substituting ( |i9|) in the asymp- 
totic evolution equation for the second order trace-free perturbations we obtain 

X$" - ~*S' + * 2 *S = Stf (*) + 0(r) , (20) 



where S a p is a function of k given in [pj| , and, as F above, is quadratic in first order 
quantities. This equation can be solved to give 

X% = ^^ + ^{kr cos kr - sin kr) + (21) 

d a p{kr sin kr + cos kr) , 

where c a p and d Ql g are functions of k, which can be expanded around r = and formally 
integrated to give 

*3<r, x) - jf (^ + U (l - ^) ) «■** . (22) 

/2 s ) 

As a result, as r — > 0, x„i tends to a constant in time which depends on the asymptotic 
values of cf)^ and d a p. These values are related to the first order initial free data by setting 
the initial conditions X a ^( T o^ k) = and x a J3 (r , k) = 0, which result in 



CoflW = fc 4 4 "^ - , M*) = — H^- ( 23 ) 

So at second order there are still four independent quantities, namely Ci,C2,a a/ 3 and b a /3, 
which correspond to the first order free initial data. Second order density perturbations can 



be computed once the metric is given [II]. 

In this way, we have shown that the second order scalar and tensor perturbations asymp- 
totically approach constant values in time that are dependent on the initial conditions. This 
demonstrates that the asymptotic behavior of the first order perturbations is stable with 
respect to the second order nonlinear perturbations. 

Now as was the case for the first order perturbations (see HH), these asymptotic constant 
values cannot be removed globally by a gauge transformation. However, this can be done 
locally: by rescaling the time coordinate, one can write the asymptotic form of (0) for the 
case of the second order perturbations as 

ds 2 = -dt 2 + e 2 v / f t A Q/3 (x, y, z)dx a dx f3 , (24) 

where A a/3 is the asymptotic form of the spatial metric 7 Q/ g ~ A a/3 +B a pT 2 , and both A a/3 and 
B a p are spatial quantities corresponding to the sum of the first and second order asymptotic 
perturbation contributions. Then as shown in ||, there exists a coordinate change (valid 
inside the event horizon of a freely falling observer at y a = 0) given by 

y a = e Ht x a 

p Ht 

where y 2 = y\ + y\ + y 2 , which allows the line element ( p4|) to be asymptotically written as 

ds 2 = -(1 - H 2 r 2 )dT 2 + (1 - tfV)" 1 ^ 2 + rW, (26) 

which represents the de Sitter solution in static coordinates 0. In this sense, as in the case 
of first order perturbations, a freely falling observer will locally see the space asymptotically 
approach de Sitter in time. 

This result is somehow unsatisfactory, as it relies on a specific frame and does not provide 
a direct insight into the fate of physical, observable gauge-invariant quantities |19| (see 



also [12, 2Tj]). Therefore, in order to shed more light on this behavior, we make use of 
curvature |18| and other invariants in order to characterize the asymptotic nature of the 
spacetime. The conformal electric and magnetic parts of the Weyl tensor can be written as 



a'~ 



E% = a 2 E« = 66* + ~0g + -Kffl + (27) 

\& a MK - & - ^) + % - 1 6 ?* . 

H a , = a 2 H° = \ 1P M 1& <s + t 1& K,s) , (28) 

where ~Q% = aQ% - ^<% B% = a 2 Rf, a% = 6f - \5p and rf*^ = e a ^/^j. It is then 
possible to show that asymptotically Q% — > 0, R% —>■ const (cf. [|12| for the expression for R% a 



perturbed to second order), Hp* — > and Ep — > const, and therefore we have #g — >■ v/A/3, 
o| -> 0, i2£* -> 0, £^ -> and if| -► 0. These in turn give 

9 2 ->• A , # 2 -»• , i? -»• 4A , (29) 

a 2 -> , £ 2 -f , iT -»• . (30) 

Now given that in the case of dust acceleration is zero and vorticity vanishes asymptotically, 
then Eqs. (p9|)-(j30|) are sufficient to prove invariantly that, locally, the nonlinearly perturbed 
spacetime approaches de Sitter. 

To conclude, we have studied the asymptotic evolution of a dust plus A spacetime with 
vanishing vorticity, by considering second order perturbations of a dust plus A FLRW uni- 
verse which is asymptotically de Sitter. We have proved the local asymptotic stability of the 
latter, using curvature and other spacetime invariants. This is a general second order gauge- 
invariant result for the case of irrotational dust with a positive cosmo logical constant, as 
guaranteed by the vanishing of ap*, R% a , Ep and Hp in the background and (asymptotically) 
at first order [JTSl rS,2D| . It is well known that vorticity decays in an expanding perfect fluid 



with equation of state p = wp, w < 2/3, and one can prove that first order scalar perturba- 
tions decay for — 1 < w < 2/3 in an asymptotically de Sitter background [Pfl . Gravitational 



wave perturbations are affected by the equation of state only through the background fL5 
and therefore their asymptotic evolution and their contribution to invariants is as described 
here. Therefore our results can in this sense be said to be more general. 

We have demonstrated the validity of the cosmic no-hair conjecture in the nonlinear 
inhomogeneous settings considered here. This, together with previous exact and perturbative 
results, also supports the picture that, starting from inhomogeneous initial conditions, the 
universe as a whole may consist of homogeneous isotropic patches that have emerged from 
inflationary phases, and others that have undergone recollapse. 

Non-linearities in the relativistic cosmology of the early and late universe have recently 
been investigated in a number of settings (see e.g. [pl,!^] and Refs. therein). We believe 



that such non-linear approximate methods are of great potential value in tackling many 
interesting problems in this field. 
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